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, Abstract. We study the dynamics of extended test bodies in flat Friedmann- 

' Robertson- Walker spacetimes. It is shown that such objects can usually alter their 

inertial mass, spin, and center-of-mass trajectory purely through the use of internal 
deformations. Many of these effects do not have Newtonian analogs, and exist despite 
' the presence of conserved momenta associated with the translational and rotational 

symmetries of the background. 

CN ■ 1. Introduction 

> ■ 

O . It is well-known that sufficiently small objects tend to fall along the geodesies of 
^ ! whichever spacetime they inhabit. While it is difficult to precisely state the limits 

i of this result (see e.g. [H El El HI E], E] ) , some qualifications are clearly necessary. Even 
O I in Newtonian gravity, the details of an extended test body's internal structure can 

i cause its center-of-mass trajectory to diverge from that of an equivalent point particle. 
^ I Inhomogeneities in the external field effectively couple to the higher multipole moments 

^ i of the mass distribution to modify the center-of-mass motion. 

^ I In principle, this effect allows piloted spacecraft to partially modify their trajectories 

simply by rearranging internal masses. A particularly elegant example of this is a 
strategy whereby artificial satellites can change their orbital parameters by cleverly 
manipulating tethered masses Different parts of the body effectively "push" or 
"pull" on local gradients in the gravitational field. In the relativistic context, test bodies 
interact with the background spacetime using their full stress-energy tensors rather than 
just their mass distributions. The control space available to alter trajectories using 
extended body effects is therefore greatly enlarged. This has a number of interesting 
consequences if very large stresses and internal momenta can be maintained. 

It allows, for example, a test body which starts at rest in a homogeneous (but 
nonstationary) spacetime to accelerate purely through the use of internal manipulations. 
No rocketry of any sort is required. Even though there is no field anisotropy for the 
body to push off of, it may still control its trajectory to some degree. This is because 
a spacetime which may appear instantaneously uniform to an observer comoving with 
the center-of-mass line needn't have this property for the various frames associated with 
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each of the body's constituents. Sufficiently large internal velocities therefore allow 
a temporal asymmetry in the geometry to have an effect very similar to a spatial one. 
Related to this is an ability for an initially nonrotating body to spin itself up by reducing 
its center-of-mass velocity (appropriately defined). 

There are other qualitatively new effects which occur in the relativistic theory as 
well. Among these is the fact that a body's rest mass is usually not conserved. This 
may be traced to changes in what may be interpreted as the energy of the system. 
The "gravitational potential energy" can vary with a body's shape or the nature of the 
applied field, for example. This affects its mass. 

We show how to make these claims precise, and illustrate them using the example 
of an uncharged extended test body inhabiting a spatially-flat Friedmann-Robertson- 
Walker (FRW) spacetime. This provides one of the ffist specific treatments of higher 
order finite-size effects in the literature. Most previous work has considered fully 
relativistic extended test bodies only as spinning point particles [211 [251 [261 [27] . These 
results assumed that quadrupole and higher order effects were negligible, although this 
is not always physically reasonable. There are many systems in which quadrupole effects 
dominate over those associated with a particle's spin. While the general forms of these 
corrections have been considered by several authors [TOl [Til [I2l [131 El IIHl [H! > th^y do 
not appear to have been previously studied in any specific systems. 

2. Laws of motion 

The constructions used here derive from an extensive formalism developed by Dixon 
[T5l [T6l [T7l [TH] and others [HI [20l [21] to generalize useful concepts in Newtonian 
gravitational mechanics to curved spacetimes. It considers the behaviour of an extended 
body described by a nonsingular stress-energy tensor T"'' with spatially-bounded support 
W. To briefly summarize, it has been shown that it is possible to convert (without any 
approximation) the laws of motion 

VaT"^^ = (1) 

into a handful of ordinary differential equations. These act on objects which may be 
interpreted as the body's linear and angular momenta. Simply put, a clear separation is 
made between those components of T"^ which are affected by stress-energy conservation, 
and those which are purely constitutive. These latter quantities may be identified as 
the body's higher multipole moments (starting at the quadrupole). A number of other 
phenomenologically interesting quantities naturally arise in this formalism as well. These 
may be interpreted as a body's inertial mass, gravitational potential energy, and so on. 

For the present purposes, Dixon's formalism is most useful when supplemented with 
a center-of-mass definition. Under mild assumptions, it is possible to define a worldline 
r which may be shown to be uniquely defined, timelike, and within the convex hull of W 
[20] . It also reduces to the standard center-of-mass definition in the appropriate limits. 
Rather than solving for the full stress-energy tensor throughout the worldtube, we only 
track F and a few supplementary quantities necessary to compute it. 
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Before proceeding with this, it is first necessary to define a body's hnear and angular 
momenta. Taking cues from the standard flat-spacetime definitions of these quantities 
[221 [23] . there should be some sense in which the momenta depend on time. At any 
particular instant, it is reasonable to define them by integrating the body's stress-energy 
tensor over an appropriate spacelike hypersurface. To be more specific, suppose that 
the center-of-mass worldline is parameterized by some function 7(s). s then defines a 
natural notion of time. A spacelike hypersurface which contains 7(3) may now be 
associated with each point on T. The set of all such hypersurfaces will be assumed to 
foliate the body's worldtube. 

With these constructions in place, note that a particle's linear momentum should 
be a vector field. In particular, let be a vector at 7(5). Similarly, the angular 

momentum S"'^{s) is taken to be a rank-2 skew-symmetric tensor at 7(s). There are an 
infinite number of reasonable definitions for momenta with these properties. A unique 
choice may be made by supposing that any Killing vector which a spacetime may possess 
implies that some linear combination of p"' and 5"* is conserved. It is clear that any 
Killing vector C,"' will imply the existence of a conserved quantity 

= / dE„(6T"'') . (2) 
We now suppose that C(^) may be constructed from the momenta according to 

Cio = P''is)Ulis)) + ls^\s)VaUl{s)) . (3) 

This uniquely defines p°'{s) and S"'^{s) in terms of T"''(x), 7(5), and The resulting 

expressions are completely independent of so they may be adopted even in the 
absence of any symmetries. Their detailed forms are not needed here, and may be 
found in UnillH]. 

Temporarily assuming that there again exists a Killing vector ([3]) indicates that 
some linear combination of the force F"(s) and torque N"'''{s) = A^[°*](s) should vanish. 
Differentiating ([3]) with respect to s, one easily finds that 

F'^^a + ^N'^'Va^b = , (4) 

where 

F'^ := 6pyds - ^^S-^i?,,/ (5) 

j^ab ._ ss^^/ds - 2p["t;^] , (6) 

and = 7"(s) = Sj'^/ds. Note that the force automatically excludes the "pole- 

dipole" component of 6p"'/ds, while the definition of the torque removes the term 
responsible for Thomas precession. Both ([5]) and ([6]) make sense even in the absence of 
any Killing vectors, so we consider these definitions to be general. 

Expressions for the force and torque may be derived from in terms of T"^ [T7l[T8] . 
If the spacetime geometry inside each slice ly nS(s) is sufficiently smooth (as discussed 
more precisely in [I5]), the resulting equations may be expanded as asymptotic series 
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involving successively higher multipole moments of the stress-energy tensor. To lowest 
non- vanishing (i.e. quadrupole) order, it may be shown that [16], [T71 118] 

Pais) ~ -lj'""'^is)\/aRbcdAl{s)) (7) 

N-\s) ^ ^r'^n-(s)R'^jUiis)) . (8) 

jabcd _ j[af)][cd] _ jcdab jg quadrupole moment of the stress-energy tensor. It is 
defined precisely in [17], although we simply note here that it does not satisfy any 
differential equations as a consequence of ([1]). All of the higher multipole moments are 
constructed such that they share this property. Their evolution equations are almost 
completely free, and depend only on the type of matter under consideration. 

Up until now, we have been assuming that 7(5) and its associated hypersurfaces 
are known. These will now be chosen uniquely using "center-of-mass" conditions. First 
let n°-{s) be a unit timelike future-directed vector field defined on F. In a sense, the 
hypersurfaces are chosen to essentially be hyperplanes orthogonal to n"-{s). Their 
precise definition may be found in [T6l [TTj [T8] . Regardless, our first center-of-mass 
condition demands that the linear momentum be proportional to this new vector field: 

= mn" . (9) 

m(s) > is some (not necessarily constant) proportionality factor which we interpret as 
the object's inertial mass. This relation suggests that be called the body's "dynamical 
velocity." This is generally distinct from the "kinematical velocity" v"". 
Next, suppose that 

Pa^"' = . (10) 

This reduces to the standard center-of-mass condition in flat spacetime. It also leaves 
the angular momentum tensor with only three independent components. Indeed, S°'^ 
may now be written in terms of a single vector S"" satisfying UaS"" = 0: 

gab ^ ^^bcd^^g^ ^ (^^^) 

Under mild assumptions, ([9]) and ( ITOi) uniquely determine n" and F throughout any 
given object [20]. We call the worldline obtained with this method the center-of-mass 
line. 

Writing down equations of motion in the standard way now requires an evolution 
equation for 7(5). Despite the highly implicit nature of the center-of-mass relations, it 
is possible to derive an exact expression for v"" — n"". Before describing this, it will first 
be convenient to fix the scale of s such that v^'Ua = —1. In general, the time parameter 
used here therefore fails to be the proper time of an observer moving on F. It is often 
very close, however. 

Defining the left and right duals of the Riemann tensor as 

^abcd • 2 cd-l^abpq i ^abcd ■ 2 ^pqcd i \^^) 
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it can now be shown that 
(m^ + *Rl^^n'S'n''S^) (n*^ - t;" + r"^) = e'^'^^'^UpSg [Ft - R^^j {n' + t") n^S^] . (13) 
is that component of the torque which appears in the expansion 

j^ab ^ ^^bcd^^j^^ + 2mr["n^l . (14) 

This decomposition is unique if UaN"- = UaT"" = 0. 

The quadrupole moment J"-^'^'^ can also be written in terms of simpler tensor fields. 
Let [191 



jabcd ^ gabcd _ ^[a^b]cd _ ^[c^d\ab _ ^^[aQb][c^d\ _ ^^^-^ 

Q'^^{s) = Q''"^^(s) may be interpreted as the quadrupole moment of the mass distribution 
seen by an observer at 7(s) with velocity n°(s). Similarly, vr"*'^ = vr"'^^! and 5''^*^'^ = 
g[ab][c(i\ _ ^cdab esscutially the body's momentum and stress quadrupoles. These 
objects satisfy 

riaQ"'' = UbTT^'"' = UaS^^"^ = , (16) 

and 

vrt'^^^l = . (17) 

There are therefore 8 independent components of tt"'^'^. Q"-^ and 5'"'"^°' each contain 6 
components, leaving 20 for the full quadrupole moment of the stress-energy tensor. This 
total could also have been deduced immediately by noting that J"*'^'^ shares the same 
algebraic symmetries as the Riemann tensor. 

The equations discussed in this section form the complete laws of motion given 
by Dixon in the quadrupole approximation. They are not the only such equations in 
the literature, however. Various alternatives have been proposed [101 [HI [El [131 [T^ - 
although each of these has drawbacks described in the introduction of [17]. To 
summarize, some contain calculational errors or inconsistent approximations. Others 
do not include the "full" quadrupole moments of the stress-energy tensor. Rather, 
some notions of the momentum and stress quadrupoles are (explicitly or implicitly) 
neglected. Some also imply undesirable results even in flat spacetime: e.g. changing 
masses or accelerating center-of-mass worldlines. These issues were all resolved by 
Dixon's formalism. For this reason and others, it will therefore be adopted universally 
for the remainder of this article. 



3. Conservation laws 



The previous section has developed enough of Dixon's formalism to allow the study of 
test body motion in the quadrupole approximation. The spatially-fiat FRW spacetimes 
we consider admit at least six linearly independent Killing vectors. Each of these implies 
the existence of its own conserved quantity, which leaves four degrees of freedom which 
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are free to affect the evolution of T. Suppose that we choose coordinates such that the 
metric takes the form 

ds^ = -dt^ + a^{t){dx^ + dy^ + dz^) . (18) 

The scale factor a{t) will be considered given. If Einstein's equation is assumed to 
hold, it depends on the averaged density and pressure of the background matter in the 
standard way. 

Generators of the isometrics associated with this FRW background may be written 
down by inspection. First define the orthonormal tetrad frame 

e(o) = dt , e(i) = a-^^di , (19) 

where i = 1,2,3. Using a 3-dimensional "vector" notation to denote triads, the minimal 
set of Killing vector fields associated with (JT8i) may now be written in the form 

f=ae, u = RxC (20) 

Here, R has been defined by the triplet {x,y,z). The cross product notation used 
to write down the rotational symmetries u then has the standard interpretation. For 
example, rotations about the z-axis are generated by the vector field uj(^z) = xdy — yd^- 

If ( ITSl) is assumed to hold throughout W , each of these six Killing vectors implies 
the existence of a conserved quantity of the form ([2]). Neglecting the body's self- 
field essentially defines the test body approximation as it is used here, although this 
assumption is not trivial. At first glance, it would seem to imply that the cosmological 
fluid penetrates the body, and that is negligible compared to the stress-energy 
tensor of the background. Most interesting systems do not satisfy these constraints. 
Intuitively, though, they shouldn't have to. This problem is a standard one, and we 
shall not attempt to address it directly. It will instead be assumed that the methods 
used here can be rigorously justified for a wide range of reasonable objects. 

With this caveat out of the way, the constants associated with the background 
symmetries - usually referred to as the system's conserved momenta - may be directly 
related to Dixon's analogs of these quantities. Slightly abusing the notation, the 
collection of scalars associated with ^ through ([21) will be referred to as the body's global 
linear momentum P = (C(C(^)), (^(c,^)))- In contrast, p"" may be thought of as the 

system's local linear momentum. Directly evaluating ([3]) using the translational Killing 
vectors given in (l20l) provides a direct relation between these two objects. Simplifying it 
using iQ, ffTTj) . and the various normalization and orthogonality conditions shows that 

p = a{mn-Hn x S) . (21) 

Here, n and S denote the triad components of the appropriate vector fields with respect 
to e. As is standard, the Hubble parameter H in fl2T!) is defined by 

(-) 

The system's two notions of angular momenta may be understood in a very similar 
way. Let those constants associated the rotational Killing vectors uj be denoted by L. 
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As before, a direct evaluation of immediately relates this object to p°-, S°-, and 7. 
After simplifying with fl2T|) . it reduces to 

(1 + \n\^j ' 

The dot product notation used here is just the standard Euclidean definition. The norm 
I rap is then n ■ n. In deriving fl23l) . the Killing vectors must be evaluated at the center- 
of-mass position 7(5). We have therefore replaced i? by 7 = (7'^,7^,7^). Notational 
conventions aside, the left-hand side of ( l23l) may naturally be interpreted as the "spin 
component" of the global angular momentum. In most cases, it differs very little from 
the local spin angular momentum S. Note, however, that it does not necessarily remain 
conserved. P and L are always fixed. 7 is not. 



4. Forces and torques 

We now have explicit relations between the global and local notions of momentum. 
Using the background tetrad, they may be used to write and S*" in terms of 7, m, 
H, P, and L. This does not uniquely fix a body's center-of-mass position, however. 
Its velocity v"' satisfies f|T3l) . which also requires knowledge of the force and torque. In 
the FRW spacetimes considered here, A^"^ is entirely determined by the vector field r" 
defined in (fT^ . With the exception of the temporal component F^-Vat, this quantity 
also fixes the force. 

These statements follow from considering the consequences of the conservation laws 
implied by the spacetime's symmetries. As explained in Sec. [2], each of these implies 
that a particular linear combination of the force and torque must vanish. In this case, 
there are therefore six relations between the the various components of F° and A^"'' 
which may be written down without any detailed knowledge of the system. Directly 
applying (jl]) for each of the Killing vectors in (|20|) shows that they have the form 

^ mHf , , 

(24) 



n 



2 



and 



Af = ^^5E^. (25) 



\n 



There remain four free parameters here: r and F°''W gt- These depend on the body's 
internal structure, and are essentially free. There exists some extended body which may 
be associated with almost any imaginable pair. The only restrictions are imposed by 
energy conditions and the consistency of the test body approximation. 

In general, specific relations between a body's internal structure and the force and 
torque acting on it may be found using the exact integral expressions in pTj [T8] . If the 
body under consideration is much smaller than a Hubble length, its force and torque 
may instead be approximated using multipole expansions. Although the conservation 
laws used to derive ( |2^ and ( l25i) are valid in the exact theory, they remain exactly 
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satisfied in any multipole approximation [12] . It is therefore fully consistent to consider 
the global momenta fixed while using the quadrupole expressions for the force and torque 
given in Sect. O 



5. Converting linear to rotational motion 

The two conservation laws (pT]) and (123|) may be used to write n and S in terms of P 
and L — ^ X P. The first of these quantities is always conserved, while the second is not. 
7 X P may vary as the body moves, which implies that the "global spin" may change. 
This suggests that moving bodies are able to change their local spin S"^ by appropriately 
manipulating their internal structures. 

Roughly speaking, the dynamics of the background spacetime mixes the momenta. 
P depends on both S and n, for example. This implies that if n or p = mn can be 
changed, so will S. Taking the norm of fl2Tl) illustrates this directly: 

Ipl = , . (26) 



1 + \HS±/m\'^ 

Here, S± represents the component of S orthogonal to n. A larger spin may therefore 
be obtained by decreasing an object's (local) linear momentum. 

This equation only takes into account the conservation of linear momentum, so one 
might expect fl23|) to rule out any such effects. It does not. Suppose for simplicity that 
P ■ L = 0, but P 7^ 0. Then n ■ S = regardless of the particle's internal dynamics. It 
follows that 

S^S.^^^. (27) 

In this case, the magnitude of the angular momentum vector coincides with that of the 
triad components displayed here: S^-Sa = S ■ S. This clarifies the interpretation of fl26|) . 

Combining fl2^ with fl7r|) allows either l^l or \n\ to be expressed in terms of m, H, 
\P/a\, and |L — 7 x P|. Although straightforward to derive, the resulting expressions 
are quite lengthy. We therefore note only that to lowest order in P, they show that 



l^^l = |L-7 X P| 



l+o{\P/ma\^) . (2^ 



This can also be deduced directly by noting that \n\ ^ \P/ma\ in this limit. The 
denominator of (127|) may therefore be expanded directly in the global momentum. (l28ll 
follows trivially. 

Regardless, this illustrates that the magnitude of a given particle's spin depends on 
its center-of-mass position. In the simplest cases, d7/ds will be nearly proportional to 
P, implying that \S±\ remains constant. More generally, though, extended-body effects 
allow significant variations in the spin. 

It is interesting to ask whether these changes are generic. Consider a particle which 
does not rotate, so 5* = 0. In this case, fl2T|) and fl23|) imply that 

7 X P = L , n = P/ma . (29) 
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Differentiating the first of tliese expressions witli respect to s sliows tliat vxn = 0. Tliis 
condition is botfi necessary and sufficient to ensure tliat tlie local angular momentum 
vanishes, v must therefore remain proportional to n. Using ( |T3|) finally shows that 
T (X P. As might have been expected directly from fl25|l . a nonspinning particle can 
only remain in that state if its internal structure is arranged such that the temporal 
torque remains proportional to the conserved momentum. Other possibilities alter the 
body's direction of travel, and inevitably impart a spin in the process. 

This statement assumes that S initially vanishes. Cases where t oc n always have 
special significance, however (at least when P ■ L = 0). Differentiating (1271) directly, it 
is clear that dS^/ds involves terms proportional to both S± and v x P. The second of 
these is more interesting. 

The center-of-mass velocity may be found using ( fT3l) and ( !24l) . Here, choosing 
T oc n implies that 

n^-v'^ + T'' (X -e^''P''npSqRlafn''n'^S^ . (30) 

It may then be shown that v (x n. Computing the spin therefore requires knowledge of 
n X P. This is most easily obtained by using fl2T|) and fl23l) to show that 

^ P + {H/m)PxS^ 

amn = . (31) 

l + \HS^/m\^ 

It follows that whenever f (xn, 

dS^/ds oc . (32) 

The direction (but not necessarily the magnitude) of the local angular momentum 
therefore remains fixed in all such cases. 



6. Zero- momentum particles 

The detailed behaviour of extended bodies with arbitrary momenta can be quite 
complicated. Suppose for simplicity that 

P = . (33) 

This is the unique condition which locks a particle's motion to that of the background 
fluid in the monopole and dipole approximations. Interesting effects appear at the 
quadrupole level, and are surprisingly rich even in this simple case. 

The spin behaviour discussed in the previous section is not retained, however. 
Dotting (|2T]) with n implies that p must vanish whenever P does. It then follows 
from (l23l) that 5* = L in this case. The local and global angular momenta agree exactly, 
so the local spin is conserved. The choice ( l33l) therefore isolates those extended body 
effects which are as independent of the angular momentum as possible. 

Contracting ( fT3l) with e^Q-j and recalling that UaT"- = now shows that up to 
a possible additive constant, the worldline parameter s is simply the cosmological 
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comoving time t. The remainder of (IT^ then yields the equation-of- motion 



d7 
m- 



rriT + H 



LxT+ {H/m){L-T)L 



(34) 



dt a(l + \HL/m\'^) 

This expression reduces to v = ad-^ / dt = f when L oc t. In other words, the coordinate 
velocity matches the coordinate (not frame) components of t"-. Unless the angular 
momentum is extremely large, this will remain an excellent approximation for the center- 
of-mass velocity of any system satisfying P = 0. 

The temporal torque r - and therefore the velocity v - may be controlled by a body 
allowed to manipulate its internal structure. Combining ([T^ . f[T^ . and ([HD, it can be 
shown that in the quadrupole approximation, 

mf ~ '^HU . (35) 

Here, H := dH/dt = dH/ds and 

H" := 7r^% . (36) 

This trace of the momentum quadrupole may point in any direction whatsoever 
(orthogonal to ra"), so an object with vanishing global momentum could have complete 
control over its direction of motion. The magnitude of this quantity limits the magnitude 
of the center-of-mass velocity. Indeed, 

1/2 

(37) 



1 + {H/mf{L 



2(T 



1 + \HL/m\'^ 

It follows that \v\ < |r| in the quadrupole approximation. These two quantities coincide 
only when S'xn = OorS'-n = 0. In other cases, a nonzero spin will reduce the 
magnitude of the particle's "drift" velocity. 

f depends on the mass parameter m, so v does as well. In general, ([5]) and may 
be used to show that 

(1 777 

— = K.^n^'n'^S" - UaF^ , (38) 
Adapted to the present case, this reduces to 

Mass changes therefore depend only on the quadrupole moments of the body's mass 
and stress distributions. They are independent of the momentum quadrupole. This 
quantity only affects the particle's velocity with respect to the cosmological rest frame, 
which is itself independent of Q"''' and 5'"*'^'^. Each of these three components of the 
full quadrupole moment may be varied independently. By construction, stress-energy 
conservation does not require any coupling between them. In principle, a properly- 
engineered spacecraft could therefore control its mass and velocity separately. 

The coefficients in front of the moments in (139|) are written so as to make this 
equation trivial to integrate if the traces Q'^a and S°'''ab remain constant. That this is 
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possible is not an accident. In the exact theory, the mass may be written as a sum of 
terms which are naturally interpreted as the total internal energy mint, the gravitational 
potential energy and the rotational kinetic energy Ej-ot of the body [18]: 

= "^int + $ + ^rot ■ (40) 

The rotational contribution here may be defined in terms of S"' and an angular 
velocity derived from an inertia tensor based on Q"'' [161 [18]. If the quadrupole moment 
remains fixed in an appropriate sense, the fact that SaS"" = cannot change ensures 
that Ej-ot remains constant. 

More generally, it is probably most straightforward to study the behaviour of the 
sum mint + Ej-ot = m — $. This may be found by noting that in the quadrupole 
approximation, the gravitational potential energy reduces to [IE] 

$~$0 + ij'^''='^i?abed, (41) 


where $o is an arbitrary constant. Applying this to the present case, one sees that 

$ ~ $0 - i(if + H^)Q\ + \h^S''\, . (42) 

Comparison with fl39l) now shows that the internal and rotational energies must evolve 
according to 



d(m-$) 1, . , ^,,dQ\ l^2dS' 



ah 



dt 2^ ^ dt 3 dt ^ ^ 

This effectively describes the rate at which energy is inductively absorbed from the 
gravitational field. 

It has been proven in the exact theory that mint remains constant for any "rigid" 
body; i.e. one whose multipole moments do not change with respect to an appropriate 
co-rotating tetrad [TS]. In this special case, full rigidity is not required even of the 
quadrupole moment. The internal energy is conserved for any object where Q"'a and 
S"'^ab remain constant. This condition is weaker than one requiring that all tetrad 
components of J'^'^'^'^ remain constant. 

When evaluating (137]) and (!39|) in specific spacetimes, it is useful write H and H 
in terms of quantities with more direct physical interpretations. For the metric fll8p . 
Einstein's equation is equivalent to [23] 

3H^ =87rp + A (44) 
3a/a = -47r(p + 3p) + A . (45) 

p and p are respectively the density and pressure of the background matter, while A 
represents the cosmological constant. Suppose for simplicity that the cosmological fluid 
obeys an equation of state of the form 

p = wp . (46) 

Also introduce the density parameter 

n := 8ttp/3H^ = 1 - A/SH"^ . (47) 
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Ignoring the particle's spin, (157|) and fl5^ now reduce to 



(48) 



and 



m ~ — 



(49) 



If Q vanishes or w = —1, the spacetime is de Sitter. In these cases, both an object's mass 
and position remain fixed for all time regardless of its internal structure. This could 
have been deduced directly using the fact that de Sitter spacetimes possess an extra four 
Killing vectors beyond the six in fl20|) . Together, these ensure that all components of 
the force and torque must vanish. There exist ten conserved quantities, among which is 
m. The center-of-mass line is always a geodesic in this case, and the angular momentum 
will be parallel-propagated along that geodesic. 

Generically, though, the object's mass will almost always change. Nearly any 
imaginable object has a nonzero mass quadrupole. Standard energy conditions should 
then imply that Q°'a 7^ 0. Excluding the de Sitter case, the first term in ( l49l) can 
only vanish if w = —1/3. This possibility only marginally satisfies the strong energy 
condition, so it is unlikely to be physically relevant. 

Still, these effects are exceedingly small even in the most favorably-chosen systems. 
If the characteristic diameter of W in the center-of-mass frame is of order D{s), the 
maximum \v\ which the particle may attain by the method we've described is of order 
(D/10^° yr)^ in the present universe. This value can only become significant for objects 
on the scale of galactic superclusters. It is a best-case estimate, however. The relatively 
low internal velocities of most types of astrophysical matter makes it very unlikely that 
there exists any observable drift velocity even in these systems. 

The situation improves slightly if we no longer assume that P = 0. For 
small momenta, v can be shown to depend on the mass quadrupole as well as H. 
This eliminates the need for a complicated internal velocity distribution. All that's 
required is an appropriate mass distribution together with some bulk motion. In 
either case, though, the maximum change in the center-of-mass velocity still scales 
as (characteristic speed) x (D/10^° yr)^- The ultra-relativistic limit |P|/ma ^ 1 is 
more interesting. If the magnitude of J"-^'^'^ is kept fixed with respect to a center-of- 
mass frame, boosting the quadrupole moment appropriately leads to a considerable 
amplification of the effects we've discussed. Of course, this case isn't particularly 
relevant to astrophysical systems either. 

Mass shifts are potentially more interesting. Although the fractional rate of mass 
change will be very small, no special structure is required to maintain it (and its sign). 
For very large objects, the total change in m over a significant fraction of the universe's 
history could be significant. Suppose that Q°'a remains constant, and ignore the stress 
quadrupole. (l39l) and (jSj) then imply that in a matter-dominated universe, an object's 
mass will change by an amount 
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Given that this scales as the cube of the redshift factor between the initial and final 
states, relatively large changes may exist even though Q°'aHf ^ m. For example, a 
galactic supercluster {D ~ 10* yr today) formed at a redshift of ~ 10 might have lost as 
much as 5% of its mass by the present time. Of course, this model is very crude. Among 
other problems, there is no reason to expect that Q°'a would remain even approximately 
constant for so long. It is also not clear how the notion of mass used here - which is 
essentially inertial - relates to the gravitational masses so fundamental to astronomical 
observations. At least in stationary spacetimes, these two types of mass are not identical 
[2T] . Understanding the analogous relations in the present case would require further 
investigation. 

Mass loss effects have been noted before in cosmological spacetimes for freely-falling 
point particles endowed with a scalar charge P |9]. In that case, the effect is due to an 
object's own scalar radiation reflecting back towards itself. This arises from the failure 
of Huygen's principle in curved spacetimes. Such exotic mechanisms are apparently not 
required to induce mass changes in freely-falling particles (even rigid ones). 

7. Discussion 

Despite the simplicity of the systems studied here, a number of interesting extended- 
body effects were found to exist. Even in the presence of linear and angular momentum 
conservation laws, it was shown that bodies can control the magnitudes of their center- 
of-mass acceleration and spin using purely internal processes. This illustrates some of 
the richness of finite-size effects in general relativity. Similar phenomena surely exist 
in other spacetimes, and their magnitudes are likely to be much larger for objects with 
realistic dimensions. 

It's also worthwhile to note that the concept of controllable motion emphasized 
here doesn't occur in simpler approximation schemes. If the quadrupole and higher 
multipole moments of a body are ignored in its laws of motion, one recovers the 
Papapetrou equations (occasionally called the Mathisson-Papapetrou or Mathisson- 
Papapetrou-Dixon equations). These follow from ([T]), and involve no free parameters 
once a definition is fixed for the center-of-mass. The motion of a spinning test body in 
a given background is completely determined by its initial conditions. 

This no longer remains true once the effects of higher multipole moments are 
taken into account. Dixon has defined these moments in such a way that stress- 
energy conservation does not affect their evolution. Their time-dependence is made 
unique by specifying constitutive relations which depend on the type of material under 
consideration. These are essentially equations of state, and are simplest to consider in 
the case of an ideal programmed mechanical device. The internal structure may then 
be considered a given function of time. 

Applying the equations of motion we've derived to more standard "passive" 
materials can be relatively complicated. This requires specifying an evolution equation 
for the quadrupole moment in terms of the various other bulk quantities relevant to 



Extended-body effects in cosmological spacetimes 



14 



the motion. For most realistic materials, all of the multipole moments couple to each 
other. It may not always be possible to ignore these couplings, in which case the 
body is probably best modeled using the full partial differential equations of continuum 
mechanics. 
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